Odd-Even Staggering of Nuclear Masses: Pairing or Shape Effect? 
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The odd-even staggering of nuclear masses was recognized in the early days of nuclear physics. 
Recently, a similar effect was discovered in other finite fermion systems, such as ultrasmall metallic 
grains and metal clusters. It is believed that the staggering in nuclei and grains is primarily due 
to pairing correlations (superconductivity), while in clusters it is caused by the Jahn- Teller effect. 
We find that, for light and medium-mass nuclei, the staggering has two components. The first one 
originates from pairing while the second, comparable in magnitude, has its roots in the deformed 
mean field. 
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The odd-even staggering (OES) of binding energies 
has been observed in several finite many-fermion systems 
such as nuclei Q], ultrasmall superconducting grains 
and metal clusters || . It manifests itself in the fact that 
the binding energy of a system with an odd particle num- 
ber is lower than the arithmetic mean of the energies of 
the two neighboring even-particle-number systems. 

In atomic nuclei, the OES is usually attributed to the 
existence of nucleonic pairing correlations . A similar 
scenario has been proposed for metallic grains (see Refs. 
HQ and references quoted therein). In both cases, the 
finite-size effects are important, and the Cooper pairing is 
well described in terms of the parity-number-conserving 
quasiparticle approach. Although the motion of electrons 
in metals is very different from that of nucleons in nuclei, 
the mechanism behind electronic and nucleonic super- 
conductivity (presence of attractive residual interaction 
which gives rise to a correlated many-fermion system) is 
indeed very similar 

So far, no evidence has been found for superconductiv- 
ity in alkali metal clusters |l0| , and the OES of binding 
energies in such systems is attributed to a very differ- 
ent mechanism. Namely, it is believed to have its origin 
in the Jahn- Teller effect which, by breaking the spherical 
symmetry of the mean field, gives rise to deformed single- 
particle orbitals Recently, Hakkinen et at [ p"3| , 
using the density-functional theory, argued that light 
alkali-metal clusters and light N—Z nuclei have a similar- 
pattern of OES, irrespective of differences in the interac- 
tions between the fermions. Hence, they concluded that 
the OES in small nuclei appears to be a mere deformation 
effect rather than a consequence of pairing. 

The main objective of this study is to analyze the phe- 
nomenon of OES in nuclei from the microscopic perspec- 
tive. Guided by self-consistent calculations, we make an 
attempt to determine and separate the pairing and mean- 
field contributions to the experimental OES. 

In the independent quasiparticle (BCS) picture 



the gap parameter, A, can be related to the binding 
energies of three adjacent systems. Assuming that the 
binding energies of even systems and those of odd sys- 
tems are locally smooth functions of the particle number 
N, the quantity 



a( 3 >(ao 



[B(N - 1) + B(N + 1) - 2B(N)] (1) 



is often interpreted as a measure of the empirical pairing 
gap. In Eq. (jl]), ttn = (— 1)^ is the number parity and 
B(N) is the (negative) binding energy of a system with 
N particles. Indeed, assuming that the ground state of 
the odd-A^ system is a pure one-quasiparticle state, one 
has A( 3 )(A0«£ fc «A, where E k = y/(e k - A) 2 + A 2 is the 
lowest BCS quasiparticle energy, e k is the energy of the 
single-particle orbital occupied by the odd nucleon, and 
A (=dB I dN) is the Fermi energy. 

Another commonly used binding-energy relation for A 
is the four-point expression [flEf, 



[3B(N - 1) - ZB(N) - B(N - 2) 



B{N + 1)} = ^[A^{N) + A^\N- 1)], (2) 



which averages the A*- 3 ) values in adjacent even and odd 
systems. In nuclei, because there are two kinds of par- 
ticles, neutrons and protons, A is calculated along the 
isotopic or isotonic chains Q . The resulting pairing gaps 
are denoted as A v and A ff , respectively. In what follows 
we argue that different physical effects determine the be- 
havior of A( 3 )(A) for odd and even particle numbers. 
Hence, the interpretation of the average (0) in terms of 
the pairing gap can be misleading. 

To investigate the Jahn- Teller component of the OES, 
we performed the Hartree-Fock (HF) calculations with- 
out pairing for nuclei with Z=9-^28 and N—Z=— 2^6. 
The HFODD code (vl.75) |6j and two different Skyrme 
parametrizations, SIII [|l7j and SLy4 |l8| , were employed. 
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For each Skyrme force, two sets of calculations were car- 
ried out: either requiring time-reversal symmetry to be 
conserved or not. (The time-odd components appear in 
the mean fields of odd and odd-odd systems.) Since, ac- 
cording to our calculations, the time-odd terms do not 
affect qualitatively the OES, the results presented in this 
work were obtained by neglecting these terms. 
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FIG. 1. Al 4) (a) and Al 3) (b) calculated in the SIII HF 
model without pairing as a function of neutron number. 

Figure [l]a displays the theoretical values of aJ> 4 ' ob- 
tained from binding energies in the SIII model. The re- 
sults obtained with the SLy4 parametrization are strik- 
ingly similar, in spite of the fact that the total binding 
energies predicted by these two forces show large differ- 
ences, approaching 4MeV in some cases. This result sug- 
gests that, unlike the total binding energy, A< 4 > weakly 
depends on the effective interaction. 

Results shown in Fig. |l]a demonstrate that the self- 
consistent mean-field theory without pairing does indeed 
predict the OES according to the criterion (||). The ef- 
fect is sizable: theoretical values of reach 30% to 
50% of the empirical OES and they, on average, gradu- 
ally decrease with mass. A rather complicated pattern of 
can be easily explained by looking at values of A^ 
presented in Fig. |]b. Values of A„ are large for 7tat=+1 
and very small for w N =-l • hence the averages 
simply reflect the simple pattern of A!, . (The behavior 



30 N 



of AV shows a very similar pattern. Also much the same 
results were obtained with the SLy4 force.) 

(3) 

We are now in a position to trace the pattern of AJ, , 
shown in Fig. [l|b, back to properties of the deformed 
mean field. Indeed, Eq. ([j]) represents the second-order 
difference with respect to N, i.e., 



2tt w A (3) (AT) 
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where g(e)=dN/de is the single-particle level density. 
Consequently, in the absence of OES due to pairing, A^ 3 ' 
represents the variation of the Fermi energy with particle 
number. In the case of a degenerate shell, the Fermi en- 
ergy does not change with N (A lies on the last occupied 
level) and A( 3 )=0. The change in A takes place when 
the valence shell e n (N=2n) is completely filled and the 
higher shell, e n +i, needs to be occupied. In this case, cor- 
responding to 7Tjv=+1, d\/dNme n +i-e n , and A^ 3 ^«(e n - 
e„ + i)/2. That is, in the absence of pairing correlations, 
A' 3 ' becomes a measure of a gap in the single-particle 
spectrum. This single-particle mechanism behind the 
OES was early noticed in Ref. || and subsequently em- 
ployed in Refs. Jfl],[f^] to explain the OES in metal clus- 
ters. The alternating behavior of A^ 3 ) in Fig. [j]b comes 
from the twofold Kramers degeneracy of the deformed 
single-particle energy levels. Indeed, the spherical sym- 
metry of the mean-field potential, which gives rise to a 
(2j+f)-fold degeneracy of single-particle levels, is pre- 
served only for doubly magic nuclei. For open-shell sys- 
tems, spherical symmetry is spontaneously broken by the 
Jahn- Teller mechanism, and the ground state is charac- 
terized by the deformed mean field, cf. Refs. |l9| , |2"o| ] . 

Results of the self-consistent calculations for nuclei 
and metal clusters can be understood in terms of the 
macroscopic-microscopic shell-correction method. In this 
approach, which can be viewed as an approximation to 
the HF treatment pi] ], the total binding energy can be 

Smacro, where 



written as B = E sp — E sp 



k=l 



(4) 



is the shell-model energy (sum of single-particle energies 
of occupied states), E sp is the Strutinsky- averaged shell- 
model energy, and E macTO stands for the macroscopic 
liquid-drop energy. 

The liquid-drop contributions to the second difference 
(^) differ for nuclei and clusters. For nuclei, the main 
effect comes from the symmetry energy term ^2|. As- 
suming the symmetry energy coefficient a/=38MeV (the 
value appropriate for light nuclei p3| ) the nuclear liquid- 
drop curvature contribution to (J3|) is 38 /A MeV, while the 
contribution of the surface-energy term is much smaller. 
On the contrary, in alkali- metal clusters, the leading con- 
tribution from the liquid-drop term to ([|) comes from the 
surface energy and is negligible. Indeed, taking the typ- 
ical value of the surface-tension coefficient for the bulk 



2 



Na, <7«0.01eV/A 2 , one obtains a very small correction 
2 E macm /dN 2 a -0.15/iV 4 / 3 eV. 

In order to evaluate the curvature contribution from 
E sp , one needs to estimate the average single-particle 
level density at the Fermi energy, g(X) [see Eq. (0)]. In 
the nuclear case g (A) = 3a /tt 2 , where a is the level den- 
sity parameter @j. Experimentally, a~A/8MeV for light 
nuclei, and this agrees well with the estimate based on 



realistic potentials |24|. Since in our HF calculations the 
effective mass is less than one, the level density param- 
eter should be additionally multiplied by the effective- 
mass factor: m*/m=0.76 and 0.7 for SIII and SLy4, re- 
spectively. Consequently, the corresponding curvature 
contribution becomes — l/g(A)«— 36/AMeV, and nearly 
cancels out the liquid-drop contribution. One can, there- 
fore, conclude that the leading contribution to the HF 
values of A' 3 ) shown in Fig. 0b comes from the single- 
particle sum (0), 



<5e=Ag)(iV)« ~(l + 7TN )(e n+1 -e n ), 
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(5) 



provided that one neglects small shifts in the single- 
particle energies due to variations of the mean field (e.g., 
deformation changes) when increasing the particle num- 
ber from N to N+l. (The effect of the mean-field varia- 
tions is seen in Fig. [l]b for particle numbers 13, 19, and 
27, i.e., around magic gaps where the shape transitions 
take place.) Quantity (||) represents, therefore, the ef- 
fective single-particle spacing between nucleonic energy 
levels. 

The above conclusion is not true for alkali-metal clus- 
ters. The level density for alkali-metal clusters can be 
estimated following Ref. §]: g(X) = 0.154m e r 2 N/h 2 . By 
taking ro=2.17A (the value corresponding to the den- 
sity of bulk Na at 500°K), one obtains l/g(A)«1.4/iVeV. 
(This result is fairly close to the Fermi-gas estimate of 
l/.g(A)w2.15/iVeV, assuming A=3.23eV.) Contrary to 
the nuclear case, since the liquid-drop component of 
the second difference is very small, the smoothed-energy 
term strongly contributes to A^ 3 ^. Moreover, the typical 
single-particle splitting due to cluster deformation is of 
the order of 0.2 eV ^5|, i.e., it is very similar to the value 
of l/g(X). For example, in Fig. 4 of Rcf. the calcu- 
lated OES parameter 2ttnA^ oscillates around zero - in 
contrast to the nuclear case presented in Fig. |l|b. This is 
because the smoothed-energy contribution to the OES, 
— l/g(X), introduces a negative shift of the same order as 
term (^|). 

Since in nuclei the mean-field contribution be is very 
small for ttn=— l f both pairing and mean-field compo- 
nents of OES can be extracted from binding energies by 
using the three-point filter A^ 3 ). As illustrated schemat- 
ically in Fig. ^, values of A^ 3 ) calculated at odd values of 
N can be associated with the pairing effect |2q] , 
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FIG. 2. Schematic illustration of various contributions to 
the OES. The odd-even energy difference, A _ , is decom- 
posed into the pairing part, A, and the mean-field part, 5e. 



A u (N) = A^(N = 2n + l), 



(6) 



while the differences of A^ 3 ) at adjacent even and odd 
values of N give information about the single-particle 
spectra, 

e n+1 - e n = 2 \a^{N = 2n) - A ( ^(N = 2n + 1)1 . (7) 
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FIG. 3. Experimental values of A V (N) = Al 3) {N) for 
iV-odd (filled circles), A? 5 for N -even (open circles), and A£, 
for iV-odd (open triangles). The thick gray line indicates the 
average trend, A=12/\/A Each point represents the arith- 
metic mean over several even-Z isotones. 

The neutron pairing gaps (^|) extracted from the ex- 
perimental binding energies are shown in Fig. [|. The 
expected quenching of neutron pairing at magic (or semi- 
magic) particle numbers iV=14, 28, and 50 is clearly seen. 
(Interestingly, the minimum at iV=20 is absent.) 

(3) 

The experimental values of A u at even neutron num- 
bers are systematically shifted with respect to A„ . Since 
the differences (Q) reflect the mean-field contributions to 
the OES, they sharply peak at magic numbers. In the 
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experimental values of the magic gaps are almost 
invisible. This is so because the effects of a large single- 
particle gap and quenched pairing cancel out in the aver- 
ages (||). Since the commonly used smooth dependence 
of average pairing gap on mass number, A = 12/y/A, 
was fitted to A^ 4 \ and not to A pTj], the values of A 
are overestimated, especially for light systems. As seen 
in Fig. ||, the values of A„ in the middle of the sd and pf 
shells are fairly similar, and are not consistent with the 
A~ 1/2 trend. 
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FIG. 4. Experimental and calculated effective sin- 
gle-particle spectrum extracted by means of relation (Q) for 
N-Z=A. 

Finally, we have extracted from the data the mean- 
field contributions to the OES according to Eq. (f7|). As 
an example, in Fig. ^ they are compared with the HF 
results for the N—Z=4 chains. In spite of the fact that 
no pairing correlations have been considered in calcula- 
tions, a good agreement between the experimental and 
theoretical levels is obtained. The fact that theoretical 
and experimental energy scales agree is a consequence of 
the fact that both SIII and SLy4 reproduce fairly well the 
symmetry energy. The level bunching predicted between 
the N=20 and A=28 gaps reflects small calculated equi- 
librium deformations in these l/7/2-shell nuclei, which 
nicely agrees with a similar grouping of levels seen in 
experiment. 

In summary, our analysis does not confirm the recent 
suggestion |rJ that the OES in light nuclei is a mere 
reflection of the deformed mean field. We have demon- 
strated that the OES in light atomic nuclei is strongly af- 
fected by both nucleonic pairing and the deformed mean- 
field. A method has been proposed to extract the pairing 
contribution to the OES from experimental data. The 
experimental pairing gaps show a weaker ^-dependence 
than that obtained previously. Since the fourth-order 
mass difference A^ 4 ^ is strongly affected by the mean- 
field contribution, it is not a good measure of pairing 
correlations, at least in light nuclei. 

In our discussion we have not discussed the singular 
behavior of binding energies of the N=Z nuclei, known 
as the Wigner energy, which dramatically influences the 
binding-energy relations (@) and (g) Conse- 
quently, the OES near the N=Z line has an additional 



third component originating from the neutron- proton cor- 
relations. In particular, the OES parameter discussed in 
Ref. |l3) , which is based on binding energy differences of 
even-even and odd-odd N=Z nuclei (hence it does not 
represent any pairing gap), is strongly perturbed by this 
third component p9[ . 
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